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Abstract 

The source term in a reaction-diffusion system, in general, does not involve 
explicit time dependence. A class of self-limiting growth models dealing with 
animal and tumor growth and bacterial population in a culture, on the other 
hand are described by kinetics with explicit functions of time. We analyze a 
reaction-diffusion system to study the propagation of spatial front for these 
models. 
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I. INTRODUCTION 



React ion- diffusion systems are ubiquitous in almost all branches of physics chemistry 
and biology dealing with population growth, fluid dynamics, pulse propagation in 
nerves, chemical reactions, optical and other processes. The basic equation describes the 
dynamics of a field variable n{x,t), a function of space and time in terms of a source term 
(also known as reaction term) and a diffusion term. An important early endeavor in this 
direction is the study of self-limiting growth models of which the most well-known is the 
Fisher equation which takes into account of a linear growth and a nonlinear decay. 
The model and many of its variants have found wide applications both from theoretical 
and experimental point of view @]. A notable feature of these models is that the source or 
the reaction terms do not involve any explicit time dependence. On the other hand there 



are situations [pHTO|| where the source terms contain explicit functions of time which put a 
constraint on the growth process in the long time limit. For example, the Gompertz growth 
is a model used for study of growth of animals and tumors, where the growth rate is 
proportional to the current value, but the proportionality factor decreases exponentially in 
time so that 

^ = rnexp(— at) , (la) 

where r and a are positive experimentally determined constants. Similarly another type of 
model proposed to analyze the growth of bacterial population in culture 0] is described by 

d/Ti 

— = kntexp{-pt^) . (lb) 

Again k and (3 are positive constants required to fit the experimental data. An important 
feature of these models is that unlike the logistic growth process the asymptotic value of the 
density function n depends on its initial population. 

Keeping in view of these experimental observations it is therefore worthwhile to generalize 
the specific cases in terms of an explicit function of time 4>{t) such that we write 
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^ = (2) 

where r is a constant for the growth process and may of the type (i) = 1 for 
exponential growth (ii) = exp(— at) for Gompertz growth (iii) = texp{—j3t'^) for 
bacterial growth, etc. 

The object of the present paper is to study a react ion- diffusion system with a reaction 
term describing a class of self-limiting growth processes (0). Since in many living organisms 
concentration dependent diffusivity [^|5|,p!T|-p!5[| has been found to be essential to the modeling 



of reaction-diffusion systems we investigate the interplay of this nonlinear diffusion and self- 
limiting growth process in the dynamics. We show that the model and its variant with a 



finite memory transport ||16|-P5(] admit of exact solutions. The dependence of the rate of 



spread of the wave front on various parameters is explored. 

II. THE REACTION-DIFFUSION SYSTEM 

We consider a reaction-diffusion system with a source term describing self-limiting growth 
and with a nonlinear diffusion term in the following form: 

where D is the diffusion coefficient. Our primary aim in this section is to provide an exact 
solution of Eq.(|^). To this end we first make use of the following transformation 

n{x,t) = u(x,t) exp (^r J (l){t')dt'^ (4) 

in Eq.(^) to obtain 

We now introduce the scaled time variable r as 

T = D f f{t')dt' = G{t) (say) (6a) 

•J 



where 

f{t)=eMr fmdt'] . (6b) 
Jo 

This reduces Eq.(^ to the following form 

dujx^r) _ d \ du{x,r) \ 
dr "axT^"^' ^ dx j 

with u{x, t) = u{x, (^"^(r)) = u{x, r) where time t has been expressed as an inverse function 
G^^{t) according to Eg . ( |6al]6b|) . 

Eq.(0) is the well-known Boltzmann nonlinear diffusion equation [[l|,|26|. Now subject to 
the initial condition of a unit point source at the origin, 

n{x, 0) = S{x) = u{x, 0) = u{x, 0) (8) 

we solve Eq.(^ under the following boundary conditions 

^lim^M(a;,r) = Vr > (9) 

and 

/ + 00 
u{x,T)dx = 1 Vr > . (10) 
-oo 

Next we seek the similarity solution of the nonlinear diffusion equation (^. We make 
use of the well-known similarity transformation 



u = T ^^^v{z) and z = xt (11) 

in Eq.(0) to obtain 



On integration Eq.(|l2]) yields 



g (i / dv\ dv ^ 

dz \ dz I dz 



3 (4) + = (13) 
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Since we are interested in the symmetric solutions with v'{0) = we have put the integration 
constant zero in going from Eq.(|12|) to (|13|). On further integration Eq.(p!3|) results in the 
solution 



viz) = iA^-z^)/6 \z\ < A 
= \z\> A 

where A is a constant which can be determined from the condition (|T0D to obtain 



(14a) 



A = (9/2) 



1/3 



(14b) 



Therefore the solution of Eq. (|^) in x and r is given by 



U(x, t) = — U - 



\x\ < 



^^1/3 







\x\ > 



At'/^ 



(15) 



It is interesting to note that by virtue of the relations ( |6a| - |6BD r is dependent on r and 
(f){t) which control the growth and self-limiting factors, respectively of the source term. 
This implies that the shock- wave like behaviour with propagating wave- front at x = Xf = 
Ar^l'^ as evident from the similarity solutions (|T5|) critically depends on the reaction terms. 
Specifically, the wave front propagates in the medium with speed 



dxf _ 1 /g^y/^ 



m 



fit')dt' 



1 -2/3 



(16) 



dt 3V 2 

where f{t) is given by (|6a|) and in turn depends on the functional form of 0(t). 

We now consider two specific cases to illustrate the spatial propagation patterns, 
(i) 0(t) = 1 

For a constant value of the model suggests an exponential growth. The relation (| 
in this case can then be utilized to obtain f{t) = exp(rt) so that r = {D/r)[exp{rt) — 1] 
Putting this expression for r in the solution (|TB|) we have after using Eq.(p 



n{x, t) 



A'^{{D/r) (exp(rt) - 1)} 



2/3 



X 



{6D /r)[exp{rt) — 1] exp(— rt) 
This solution clearly has a sharp wave- front at Xf = Ar^^^ which propagates at a speed 



(17) 



'^ = \a {Dr^y^' exp(rt) (exp(rt) - 1)-'/=^ (18) 

To illustrate the spatial propagation of the population n{x,t) in time we plot in Fig-1 
the spatial shock-wave like patterns for r = 1.0 and D = 1.0. It is apparent that the sharply 
peaked distribution at t = starts spreading relatively slowly with peak at x = diminishing 
with time upto a period t = 0.1 . Beyond this time the spatial growth of population becomes 
comparatively large and it diverges due to the combined effect of exponential growth and 
nonlinear diffusion. For a much lower growth rate (r = 0.001), however, the population 
spreads monotonically due to the nonlinear diffusion which overwhelms the effect of growth 
process. This is evident in Fig-1 (b). 

(ii) 0(t) = texp(-/?t2) 

With the above expression for 0(t) for bacterial self-limiting growth we obtain from 

-eg) 



fit) = exp[(-r/2/?)(exp(-/?t) - 1)] (19) 



and 



r = D exp(r/2/3) T exp[(-r/2/3) exp{-(3t')]dt' (20) 
Jo 

By defining z = {r/2(3) exp(— /3t) the above expression can be reduced to the following form 

^exp(r/2/?) /-(rmcxpC-/?*) exp(-z) , , . 

r = -D / dz (21) 

(3 J(rm Z 



SO 



The integral in ( pTD can be put into a standard form with the help of £'i-function |^ 
that r can be expressed as 

r = /^ ^^P(^/2/^) (-r/2/5) - Ei ((-r/2/3) exp(-/3t))] (22) 

The corresponding density n{x,t) and the speed of the wave front dxf/dtaXxf are given by 

_ {D exp{r/2f3)/f3f' [Ez (-r/2/?) - Ez (-(r/2/j) exp{-f3t))f' - x' 
' (QD/p) [Ei{-r/2p) - Et{{-r/2f3) exp(-/3t))] exp[(r/2/3) exp(-/3t)] ^ ' 

and 



D ^^ j ^ [i?z(-r/2/3) - Ez((-r/2/3) exp(-/5t))]^/^ (24) 

respectively. 

In Fig-2(a,b) we show the shock-wave hke spread of population by plotting n{x,t) vs 
X for several values of time for D = 1 and r = 1 . Since (3 puts a limit to the growth at 
large time the peak of n{x,t) at x = as shown in Fig-2(a) {(3 = 0.1) does not increase too 
much as compared to the earlier case considered in Fig- 1(a). It has been observed that for 
a unique value of /? > 1.0 there is a monotonic decrease in the peak population n{x,t) at 
X = 0. For smaller values of (3 (Fig-2(b)) the spread is similar to that in Fig- 1(a). In Fig-2(c) 
we exhibit the spatial front propagation for several values of growth rate r at a time t = 1.0 
keeping D = 1 and (3 = 0.01. It is apparent that with increase of r the reaction dominates 
over diffusion so that the peak population at x = increases compared to spreading. 



dt \ 



III. EFFECT OF FINITE MEMORY TRANSPORT 



We now generalize the proposed reaction-diffusion model to include the effect of finite 
memory transport. It has been observed that an animal's movement at a particular instant of 
time often depends on its motion in the immediate past. This results in a delay in population 
flux, or a memory in the diffusion coefficient. A number of attempts have been made in 



the recent literature p!6|-p5[] to analyze the delayed population growth in several models and 
related context in heat conduction and transport processes. To consider a finite memory in 
the present model we modify the nonlinear diffusion term in Eq.(|^) to the following form: 

dn{x, t) 



dt 



rn{x, t)(j){t) + 



d_ 

dx 



rt Qflix T] 

D'-f / exp[— 7(i(: — r)]n(x, r) — — ^ — dr 

J CJju 



(25) 



where 7 refers to the inverse of relaxation time. The population flux takes into account 
of the relaxation effect due to the delay of the particles in adopting a definite direction of 
propagation. Differentiating both sides of the above equation with respect to t and using it 
again we obtain 



dn ■ d 

{rcj) — 7)— — h (r0 + r(j)'y)n + 



dn 
D'jn— 
ox 



(26) 



dt"^ dt dx 

In the limit of vanishing relaxation time i.e, I/7 — s>0 Eq.(p6D reduces to Eq.(^. When 
memory effects are taken into account, the dispersal of the organisms are not mutually inde- 
pendent. Hence the correlation between the successive movement of the diffusing particles 
results in a delay in the transport. Thus Eq.(p6|) is a typical form of a delayed transport 
equation. 

We now consider a specific case 0(t) = 1. Substitution of the traveling wave form N{z) 
(= n{x, t)) with z = X + ct satisfies 

where c is the speed of the traveling wave to be determined. 

We now consider the trial solution of Eq. (pT]) of the form N{z) = N^exp^sz'') subject 
to the initial condition that at z = 0, N = Nq, where s and b are positive constants to be 
determined. Substitution of this solution in Eq.(^) yields the following relation 

-r7] expisz'') - D^Nosb[2sbz^^^~^^ + {b - l)z^^-^^] 

X exp(2s2'') = L{z) = (28) 

For L{z) = 0, for all z, the coefficients of exp(sz'') and exp(2s2;^) within the square brackets 
must vanish identically. For this the only acceptable solution for 6 is 6 = 1. We obtain 

2s'^D-fNo = (29a) 

and 

c2s2-cs(r-7)-r7 = (29b) 

From the above two equations the solution for s is given by 

c[(l/7) - (1/r)] + [c' ((I/7) - il/r)f + 4/r ((cV7) + 2DiVo)]V2 



2 [(cV7r) + {2DNo/r)] 
8 



(30) 



In the limit of instantaneous relaxation, i.e, I/7 — >0 Eq.(|30|) yields 

-1 + (1 + {DNor/^)^' 



c 

s = — 



(31) 



Furthermore the above expression in the limit of weak diffusion D —>-0 we obtain from 
Eq.(^T]) after a Taylor expansion 

r 

s = - (32) 

c 

To determine the speed of the propagation of the wave front we now rearrange the 
solution for s in ( ^O]) to obtain 

_ (r - 7) + [(r - 1? + 4(r7 - 2s'D^No)f' 

For real values of c, the quantity inside the square root must be positive, which determines 
the minimum value of c for s = r/c [ Eq.(^) ] as 



(r + 7)^ 

Eg ■ (|27|) therefore admits of an exact traveling wave like solution 

N{z) = iVoexp ''^'^ " + ^'''^'^ ' + + 2D^No)f' 



(34) 



(35) 



2(c2 + 2D7iVo) 

It is interesting to observe that the speed of the traveling wave front not only depends 
on nonlinear diffusion and growth rate but also on the initial concentration and memory. A 
comparison of the solutions in this section and in the previous one shows that ([35| ) does not 
reduces to Eq.([T7|) in the limit of vanishing relaxation time (I/7 0) although Eq. (p6D goes 
over to Eq. (^ under this condition. This is because of the fact that the nature of the partial 
differential equation changes due to the inclusion of relaxation terms and also the boundary 
conditions for the shock wave like 'diffusing solutions' (p!7| ) are different for the travelling 
wave front solution (pSj). The nature of the two solutions are thus generically different. We 
point out in passing that the dependence on initial concentration on speed as shown in 
is rather an unusual feature in reaction-diffusion system. 
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IV. CONCLUSIONS 



In this paper we have analyzed a class of reaction-diffusion systems in which the kinetic 
term describes the self-limiting growth processes of the Gompertz type and is an explicit 
function of time. We have shown that the model can be solved exactly to analyze the spatial 
front propagation problem. To make the model more realistic we have included the effect 
of finite relaxation to concentration-dependent diffusive processes. In view of the fact that 
the source terms have their direct relevance on experimental measurement on animal and 
tumor growth or bacterial culture we think that the solutions discussed in this paper will 
be pertinent in the context of reaction-diffusion systems, in general. 
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FIGURES 

FIG. 1. Evolution of spatial front in time for the model with <^{t) = 1. (a) The population 

n{x,t) is plotted against x for different times using r = 1.0 and D = 1.0. (b) Same as in Fig. (la) 
but for r = 0.001 (units arbitrary). 

FIG. 2. Evolution of spatial front in time for the model with (j){t) = texp{—pt^). (a) The 
population n{x,t) is plotted against x for different times using r = 1.0, D = 1.0 and /? = 0.1. (b) 
Same as in Fig.(2a) but for (3 = 0.01. (c) The population n{x,t) is plotted against x at t = 1.0 for 
different r using D = 1.0 and (3 = 0.01 (units arbitrary). 
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